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A study of secondary instabilities in ac-driven electrocon- 
vection of a planarly aligned nematic liquid crystal is pre- 
sented. At low frequencies one has a transition from normal 
rolls to a zig-zag pattern and at high frequencies a direct tran- 
sition from normal to abnormal rolls. The crossover defines 
a codimension-2 point. This point is also the origin of a line 
of homoclinic bifurcations, which enable the transition from 
zig-zags to abnormal rolls at lower frequencies. 
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A large degree of universality has been found in pri- 
mary instabilities leading to periodic patterns in spatially 
extended systems [Q. A full classification was achieved. 
For example, in axially anisotropic quasi-2D systems un- 
dergoing a steady supercritical bifurcation one can have 
at threshold either normal rolls (NRs), which orient ac- 
cording to the preferred direction, or oblique rolls, which 
occur in two symmetry-degenerate variants and may su- 
perpose to give rectangles [|| . 

A full classification of the secondary instabilities that 
may destabilize the primary patterns appears rather pro- 
hibitive, except in the quasi-lD case ||^. For Rayleigh- 
Benard convection in simple fluids, the prime example for 
isotropic quasi-2D pattern formation, the relevant sec- 
ondary instabilities generate the celebrated Busse balloon 
and have been analysed in detail Q. By contrast, in the 
best studied axially anisotropic system, namely ac-driven 
electroconvection (EC) in planarly aligned nematic liquid 
crystal layers, the secondary instabilities arising when the 
main control parameter (the amplitude of the voltage) is 
increased have only recently been understood starting 
from the hydrodynamic equations [||. It was found that 
in the range where one has NRs at threshold the well- 
known long-wave zig-zag (ZZ, or undulatory, or wavy) 
instability occurs only for frequencies / (the second con- 
trol parameter) below some value /ab- 

For / > /ar one should have a homogeneous instabil- 
ity. The nematic director, which near the primary thresh- 
old lies within the plane made up of the NR wave vector 
(II £') and the layer normal (|| z), rotates out of this plane 
in one of two symmetry-equivalent directions. Since the 
director is anchored at the upper and lower plate along x, 
this involves a twist deformation of the director field, and 
one may speak of a twist mode. The resulting roll pat- 



tern was termed 'abnormal rolls' (ARs), following simi- 
lar observations in EC with homeotropic anchoring 
Since homeotropic systems can be rotationally invariant 
or nearly so (in the presence of a weak planar magnetic 
field) the AR-transition can occur at or near threshold, 
allowing for a description of the scenario in terms of cou- 
pled Ginzburg-Landau equations |^]. However, one of 
the most interesting effects elucidated in Ref. |^ , namely 
the restabilization of ARs above the ZZ-instability, is 
hardly expected under these conditions Q. 

We here present a study of the predicted crossover from 
the ZZ- to the AR-transition in planar EC. Besides iden- 
tifying experimentally this new type of codimension-2 
(C2) point, we unravel the behavior in its neighborhood 
in the voltage-frequency plane experimentally and by a 
normal-form analysis involving two relevant modes. We 
show that at this C2-point a hitherto unsuspected line 
originates, which limits from above the regime where one 
may have ZZ-patterns. At this line one has an unusual 
type of homoclinic bifurcation from ZZs to ARs. A par- 
ticularly interesting feature is the hysteresis found when 
this line is crossed from above. Then ARs persist down 
to a well-defined stability limit. The normal form shows 
that this class of two-mode scenario is of general impor- 
tance, e.g. for certain types of line defects. 

We perform experiments on EC in the conduction 
regime of the liquid crystal 4-methoxybenzylidene-4'-n- 
butylaniline (MBBA) at a temperature of 15°C. We use 
a channel geometry that is only moderately extended in 
the x-direction, namely a capacitor, where a 25 x 10 mm 
ITO-coated glass plate faces another glass plate with a 
Lx = 315 /xm wide and Ly w 10000 fim stripe of ITO cre- 
ated by an etching process j|J^ . The distance d — 24 iim 
between the two plates scales the wavelength of the con- 
vection pattern, which arises when the amplitude V of 
the driving ac-voltage passes a certain critical value Vc- 
A mechanical treatment of the surfaces (rubbing) fixes 
the planar orientation of the director almost perpendic- 
ular to the y-direction, which is defined by the lateral 
boundaries of the channel. The patterns are visualized 
by the standard shadowgraph technique ||l^ and detected 
using a microscope and a video-camera. 

The left-hand side of Fig. 1 shows the NRs in the 
low frequency regime. They arise from the unstructured 
ground state via a primary instability, and destabilize at 
a higher voltage due to a secondary ZZ-instability. The 
resulting pattern is shown in the middle part. At an even 



1 



higher voltage, ARs appear as shown in the right-hand 
side of Fig. 1. They look very similar to NRs — a spe- 
cial optical technique is needed in order to measure their 
underlying twist [ pH . This image represents the first di- 
rect experimental observation of the ZZ-AR transition — 
presumably supported by the lateral boundaries — in EC 
in a planarly aligned liquid crystal. 
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FIG. 1. Convection patterns in a 315 x 1020 part 
of the channel at 3 different ac-voltages at 43 Hz. 

The procedure to obtain the five different relevant tran- 
sitions in this scenario is illustrated in Fig. 2. The mea- 
surement of the primary instability is presented in the 
lower part. Here the light intensity modulation along 
the X-axis was analyzed by means of Fourier analysis. 
The intensity of the peak corresponding to the critical 
wave number is shown as a function of the driving volt- 
age. This number is supposed to increase linearly with 
the distance from the threshold voltage Vc- The linear 
extrapolation of the data to zero determines the onset of 
convection. In our finite channel this bifurcation is im- 
perfect . We have measured the intensity for increas- 
ing (open squares) and decreasing (solid circles) voltage: 
This transition has no hysteresis. 

The angle of the ZZs with respect to the boundaries 
was measured by means of a correlation technique, at 
fixed position where this angle turned out to be maxi- 
mal. This order parameter is shown in the middle row 
of Fig. 2. It is obvious that ZZs appear at a driving fre- 
quency of 43 Hz, but not at 80 Hz. The onset of ZZs 
at 43 Hz occurs at 6.55 V, as measured by means of a 
square-root extrapolation. The transition is imperfect 
and has no hysteresis. The transition from ZZs to ARs 
was determined with a threshold criterion (dotted line). 
It has a large hysteresis. The two points measured at the 
transition from ARs to ZZs are transients, which cannot 
be avoided within the pacing of 1 minute used here. 

In order to determine the onset of ARs we have mea- 
sured the twist angle by a special optical setup [|ll| , which 
analyzes the ellipticity of the transmitted light and will 
be described in detail elsewhere The results are 

shown in the upper column of Fig. 2. The transition to 



ARs is easy to identify at the driving frequency of 80 Hz. 
An (imperfect) supercritical pitchfork bifurcation is indi- 
cated at a voltage of 7.2 V: The first direct experimental 
demonstration of the transition from NRs to ARs in EC. 

At the driving frequency of 43 Hz the situation is more 
complicated. Here the ZZ-pattern mediates between the 
regime of NRs and ARs. The ZZs are accompanied by 
a twist. Thus the hysteresis in the ZZ-pattern manifests 
itself also in a hysteresis of the measured twist angle. 
The bifurcation from NRs to ARs can now be identified 
only indirectly, because both patterns are unstable with 
respect to ZZs in the vicinity of this bifurcation point. 
We thus measure these points by an extrapolation from 
the data obtained in the regime of ARs, assuming an 
imperfect bifurcation as indicated by the solid line. 
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FIG. 2. Measurements at two frequencies. Open squares 
(solid circles) correspond to increasing (decreasing)voltage. 
We show the twist angle (top), the zig-zag angle (middle), 
and the strength of convection (bottom). The solid lines are 
fits to determine the thresholds. 

Fig. 3 is a stability diagram, which displays all five bi- 
furcation lines discussed above. The crosses correspond 
to Vc- The open squares indicate the ZZ-instability. The 
solid squares (open diamands) mark the transition from 
ZZs to ARs (ARs to ZZs) taking place at increasing (de- 
creasing) voltage. The solid circles indicate the bifurca- 
tion from NRs to ARs. In order to simplify a compar- 
ison with the theory, the same data are shown in the 
upper part of Fig. ?? with a rescaled voltage, namely 
e — V"^ /V^ ~ 1. Note the similarity of this measurement 
with Fig. 3 of Ref. Q (those calculations were done for a 
different material) . The solid lines are least-square fits of 
data points above 67 Hz to the theory presented below. 



2 



0.20 



CO 




a 




a 


0.15 


s 








C 
■p 


0.10 






o 










0.05 




7.5 




7.3 




7.1 






ID 


6.9 


be 




CIt 

-u 
1— 1 


6.7 


o 
> 


6.5 




6.3 




6.1 




• 8 ° 



40 50 60 70 80 
Frequency (Hz) 



90 



FIG. 3. Experimental stability diagram. We show the 
threshold of convection (crosses), the ZZ-instability (open 
squares), the disappearance of ZZs for increasing voltage 
(solid squares) , the destabilization of ARs for decreasing volt- 
age (open diamonds), and the AR-instability (solid circles). 

For a theoretical description of the scenario we adopt a 
phenomenological approach based on a Ginzburg-Landau 
expansion around the bifurcation from NRs to ARs in the 
vicinity of the C2 point. Let (j) describe the (real) am- 
plitude of the twist mode, which plays the role of the 
order parameter. The pre-bifurcation state sustains a 
spontaneous periodic pattern whose phase represents a 
soft (Goldstone) mode which has to be included in the 
description. We write the phase as qx + V9, such that 
qx + V9 is the local wave vector. Allowing only for vari- 
ations along y (as is adequate for the channel geometry) 
the Ginzburg-Landau equation for the pitchfork bifurca- 
tion from NRs to ARs is given by 



dt(j) = (/i 



(1) 



with positive g and K and (reduced) control parameter 
fi. The last term describes the bias when the rolls are 
(slightly) oblique (the local roll angle is arctan{dy9 / q)) . 
We allowed for a misalignment between the rubbing di- 
rection {(f> — 0) and the x-axis defined by the channel 
geometry. 7 must be chosen positive, so reorienting the 
rolls favors rotation of the director in the opposite sense 
1^. From a fit to the experiment (Fig. 2, top right) one 
finds — 0.0011. The dynamics of the phase mod- 

ulation is governed by the equation 
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dyJ, J 



D{dy9 + Pa)-iu + h^''/3)(f>, (2) 



with positive D and h (as it turns out). The first term in 
J describes ordinary phase diffusion, and the second ex- 
pression represents the coupling to the twist mode. The 



nonlinear term is essential, because v crosses zero in the 
region of interest. 

Straight-roll solutions (normal or oblique) are char- 
acterized hy (p = 4>s {— const.) and dy9 — P {= 
const.) and their stability is derived from the growth 

rate A = 
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C of modes ^ e*^*', where 

B = + 3g(l)'i + (K + D)p^, C ^[ju-Dn + 3{gD + 
/i7/3)(^2]p2 _|_ j^j^pi YoT stabihty one needs B > 
and C > 0. Thus, for NRs with 0^ = and Pa = 0, 
there is a homogeneous instability at /i = 0, leading to 
ARs and a long-wave ZZ-instability at /j, = — Szz^ 
with Szz = l/D. For negative v one first has the ZZ- 
instability and eqs. indeed describe the observed 
crossover scenario, as indicated in Fig. ??, lines AR and 
ZZ. From the expression for C one sees that a nonzero 4> 
suppresses the ZZ-instability. For negative v this effect 
leads to restabilization of ARs {(i>\ft. = i^ls) above the 
line /i = SrsV with Srs 
Fig. 4, line RS. 
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FIG. 4. Theoretical stability diagram for S^z/Sht ~ —2/5. 
Solid (dashed) lines are for Pa ^ {-yg^^^Pa = 0.0011). 

We now analyse the scenario expected when fi is varied 
for negative ly. For static solutions of (|l|,^ one has J = 
const. Eliminating dy6 leads to 



(3) 



with g' = g(l — Szz/Shb), which is integrable. Invok- 
ing the analogy of a point particle (coordinate (f), time 
y) one sees that the bounded solutions are either con- 
stant or periodic. The channel geometry requires that 
rolls remain on the average oriented along y, so that 
the average of dy9 is zero. From (||) one then has 
J = DPa- < + \h4)'^)4) > (< ... >= spatial aver- 
age). 

First we look for ZZ-solutions in the case Pa = 0, 
where (p oscillates symmetrically around zero (J — 0). 
From (js]) one gets a one-parameter family of periodic 
solutions above the ZZ-line (apart from phase shifts). 
In particular there is the homoclinic (or actually hetero- 
clinic) limit where the solution degenerates to a widely 
spaced array of domain walls separating regions where 
4> approaches the constant solutions -ii4>zz with (pzz = 
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\/ {n — Hzz)/ g' ■ The undulations observed under increase 
of the voltage are (presumably) approximated by this so- 
lution. The maximum value of the roll angle is given by 
arctan(p^^/g) with p^^ = 30^372 (m - fJ-zzV^'^if^ - fJ-hb)- 
The roll angle (and thereby the undulation) first increases 
with fj, and then decreases, reaching zero at the line 
^J' — ^'hb — ShbV. There 4>zz coincides with 4>ar and one 
is left with an array of marginally stable domain walls 
separating AR-regions with alternating director twist. 
Above this line (and also for positive v) only moving do- 
main walls exist stably (for details see jl3|). 

Therefore, when the line HB is reached, the domain 
walls annihilate pairwise, and a single-domain AR is 
established. From now on ARs persist under changes 
of the parameters until their stability limit is reached. 
The integration constant J maintains the AR value 
— {v + ^h(j)\jj)(j)AR (= at Hhb)- Thus, on lowering /i, 
ARs persist down to the line RS where a discontinuous 
ZZ-instability occurs, as observed. The twist angle in the 
ZZ-state is larger than in the coexisting ARs, in agree- 
ment with the experiment (Fig. 2, top left). The solid 
lines in Fig. 3 represent a five-parameter fit of 4 lines 
through a common intersection point (the C2 point) with 
slopes incorporating the relation Srs = 35 /sit -2 ' ^^^^ 
relation is in fact invariant under the linear mapping that 
connect the control parameters /i and i' of the model with 
the experimental control parameters e — eAR and f ~ fAR- 

The observed asymmetry of the undulations can be 
accommodated by choosing Pa ^ 0. In Fig. 4 the cor- 
responding scenario is shown (dashed lines). The ZZ- 
instability remains sharp, since translation invariance 
along y remains intact. At the homoclinic bifurcation 
again J = 0. Now the length of the shorter arms of the 
ZZs vanishes there. The selected wavelength of the ob- 
served ZZ undulations presumably results from the finite 
width of the channel [p^ . 

The nature of the hysteresis found is actually quite 
unique. One ingredient, typical for pattern-forming sys- 
tems, is the conservation of the number of periodic units 
(here ZZ-undulations) due to pinning at the boundary 
(here the short ends). Another ingredient is the unusual 
homoclinic bifurcation of ZZ patterns from ARs. We 
note that the ZZ-instability comes out naturally, which 
underscores its importance for EC. In fact, for many 
other nematics, like the one used in Ref. g, the ZZ- 
instability line joins the primary bifurcation line at an- 
other significant C2-point, which separates the regimes 
where oblique and normal rolls appear at threshold. 

In summary, two nontrivial qualitative features of the 
model are in accordance with the experimental data: A 
hysteretical transition between ZZs and ARs vanishes at 
the C2-point, and the width of this hysteresis decreases 
linearily with the distance from this point. 

The model should be applicable to other instabili- 
ties exhibiting a similar symmetry. Candidates are line 



defects and some types of domain walls in the bend- 
Freedricksz distorted state in nematics. Sometimes in 
these ID extended structures the director can escape out 
of the symmetry plane, which may mimic the transition 
from NRs to ARs described by (f>. If the position of the 
line/wall is not fixed from outside it can be described by 
our phase variable 9. In those cases where one has a po- 
tential (no dissipative driving) our model predicts ^yv < 
so that the ZZ-instability always occurs first, which ap- 
pears to be consistent with experiments [Q. The model 
is not applicable if the coupling terms between the two 
active modes vanish by symmetry, as is the case in Ising- 
Bloch-type transitions of domain walls . We also re- 
mark that the validity of the model is restricted to the 
immediate neighborhood of the C2 point. 
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